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^ i Abstract 

On 
O 

C*~) ' The successful prediction of the weak mixing angle suggests that the 

' 

effective theory beneath the grand unification scale is the minimal super- 
| symmetric standard model (MSSM) with just two Higgs doublets. If we 

c-j . further assume that the unified gauge group contains SO(10), that the two 

light Higgs doublets lie mostly in a single irreducible SO(10) representa- 
tion, and that the t, b and r masses originate in renormalizable Yukawa 
' interactions of the form I63OI63, then also the top quark mass can be 

predicted in terms of the MSSM parameters. To compute mt we present 
a precise analytic approximation to the solution of the 2-loop renormal- 
ization group equations, and study supersymmetric and GUT threshold 
corrections and the input value of the b quark mass. The large ratio of 
top to bottom quark masses derives from a large ratio, tan /3, of Higgs 
vacuum expectation values. We point out that when tan/3 is large, so 
are certain corrections to the b quark mass prediction, unless a particular 
hierarchy exists in the parameters of the model. With such a hierarchy, 
which may result from approximate symmetries, the top mass prediction 
depends only weakly on the spectrum. Our results may be applied to any 
supersymmetric SO(10) model as long as At ~ A& ~ A T at the GUT scale 
and there are no intermediate mass scales in the desert. 
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I. INTRODUCTION 



The standard model of particle physics is extraordinarily successful, describing all 
known properties of the elementary particles in terms of just 18 free parameters. Nev- 
ertheless, the model leaves so many questions unanswered that numerous ideas and 
speculations leading toward a more fundamental theory have been developed. While 
there is no hard evidence to support any of these speculations, experiments have pro- 
vided some hints. In particular, the only parameter of the standard model that has been 
successfully predicted with a high level of significance is the weak mixing angle, which 
is a prediction of supersymmetric grand unified theories (SUSY GUTs) [|J. 

While this is just a single parameter, the excellent agreement of data with the simplest 
SUSY GUT suggests that it is worthwhile to pursue other predictions of similar SUSY 
GUTs. This is harder than it sounds. The reason is that the weak mixing angle has 
a unique status within these theories: it is the only parameter which can be predicted 
by knowing only the gauge group structure of the model. In fact all one needs to know 
is that the gauge group is, or breaks to, SU(5) ||. All other predictions require 
additional, model-dependent information about the theory. A good example of this is 
the proton decay rate: it can only be computed after making an assumption about the 
spectrum of the superheavy colored states. It is a very interesting quantity, since the 
simplest possible structure for this superheavy spectrum gives a decay rate that may 
well be accessible to planned experimental searches. Nevertheless, only minor changes 
in the theory can lead to very large suppression factors in the rate. 

A potentially copious source of predictions is the flavor sector, responsible for the 
quark and lepton masses and mixings. An early success of GUTs was the prediction of 
the bottom quark to tau lepton mass ratio, mb/m T 0; however, several considerations 
make this less impressive than the weak mixing angle prediction. First and foremost is 
the low numerical significance of the m,b/m T prediction: it cannot be made with much 
accuracy as long as the strong coupling and the top quark Yukawa coupling (which 
cannot be neglected for a heavy top) are not known very well, and it cannot be compared 
accurately with experiment without better knowledge of mj. This low significance is 
especially troubling given that the theory is predicting only one of the 13 independent 
flavor parameters of the standard model. Nevertheless, with a heavy top quark, an 
acceptable value of m&/m r requires more than one light Higgs doublet — which provides 
another interesting hint pointing towards low-energy supersymmetry ||. Finally, the 
simplest flavor sector which leads to this prediction immediately fails when extended to 
the lighter generations: m s /md = m^/me is unacceptable. To overcome these objections 
it is necessary to construct a more complicated flavor sector of the grand unified theory. 
Here there is a delicate balance: more structure requires further assumptions, but these 
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are perhaps justified if there are additional predictions. This approach was developed 
long ago |pHTl[] and has received considerable attention recently Hl2| . Using the full 
power of the grand unified group SO(10) |T3j it is possible to obtain predictions for 7 
of the 13 flavor parameters. Further development of the flavor sector can also lead to 
predictions for neutrino masses pTT Jl^| . Despite these successes, one still has to admit 
that these schemes are based on the hope that the flavor sector at the grand unified 
scale is particularly simple: the quark and lepton masses must originate in just a few 
grand unified interactions. If there are many such interactions the predictions are lost. 
This is a particularly acute problem for the lighter generations. The smallness of these 
masses can be understood if they arise from higher dimensional operators. However in 
this case there is a very large number of operators that could be written down, and the 
restriction that just one or two of these operators contribute to the masses involves some 
strong assumptions. 

We are unable to completely avoid this dilemma: to obtain quark and lepton mass 
predictions from grand unified theories, assumptions about the underlying flavor struc- 
ture of the theory must be made. Of all the flavor predictions of GUTs those pertaining 
to the heaviest generation are most direct and subject to the fewest assumptions. In this 
paper we pursue a scheme for predicting the top mass which is unique in its simplicity. 
We attempt as complete and accurate analysis of this prediction as possible: the aim is 
to predict the top quark mass to within a few GeV. 

There are many approaches in the literature which result in predictions for the top 
mass. We would like to emphasize the differences between two such approaches and 
a third one which we shall take. The first is the infrared fixed point behaviour of the 
renormalization group (RG) equation for the top quark Yukawa coupling [|TJ]. This is an 
argument that certain values for the top mass are more probable than others if all GUT 
scale Yukawa couplings are equally probable. The second framework for predicting the 
top quark mass is that of textures for the generation structure of the Yukawa coupling 
matrices. The top mass is given in terms of lighter quark masses and entries of the 
Kobayashi-Maskawa mixing matrix |10|JTT[] . This necessarily requires assumptions about 
the masses of the lighter generations. Finally, in the context of grand unified theories, 
it is sometimes possible to predict the top quark mass from a consideration only of the 
heaviest generation, with the b and r masses as inputs from experiment (T^JT?]]. Consider 
the Yukawa interactions of a supersymmetric grand unified theory which lead to masses 
for the heaviest generation. There are three observable masses: m t , m?, and m T . Since the 
Higgs doublet which leads to the top mass is forced by supersymmetry to be different 
from the one which gives mass to b and r, these three masses necessarily depend on 
the parameter tan (3, the ratio of the vacuum expectation values (VEVs) of these two 
doublets. If the grand unified theory has two or more independent Yukawa parameters 
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contributing to the heaviest generation masses, then there will be no prediction when the 
heavy generation is considered in isolation. The only possibility for a prediction resulting 
from consideration of the heaviest generation alone is that the t, b and r masses originate 
predominantly from a single Yukawa interaction. This immediately excludes the grand 



unified gauge group SU(5) from consideration [ 18[ . The simplest supersymmetric [19 
grand unified gauge group which allows relations between masses in the up and down 
sectors is SO(10). 

In this paper we study the top quark mass prediction which results from the following 
three assumptions: 

(I) The masses of the third generation, m t , and m T , originate from renormalizable 
Yukawa couplings of the form I63OI63 in a supersymmetric GUT with a gauge 
group containing (the conventional) SO (10). 

(II) The evolution of the gauge and Yukawa couplings in the effective theory beneath 
the SO(10) breaking scale is described by the RG equations of the minimal super- 
symmetric standard model (MSSM). 

(Ill) The two Higgs doublets lie predominantly in a single irreducible multiplet of 
SO(10). 

We find it highly significant that such simple and mild assumptions are sufficient for 
predicting the top quark mass with an accuracy of a few GeV, in terms of a few pa- 
rameters of the MSSM (which could be measured experimentally). In Sees. II VII we 
in fact assume that the two light Higgs doublets lie completely in a single irreducible 
multiplet, while in Sec. IX we return to the effects of mixing with other multiplets. We 
will very rapidly be led to the result that the SO (10) multiplet containing the Higgs 
doublets is (almost) necessarily the 10 H , so that the relevant Yukawa interaction is 
16 3 10 H 16 3 (where 16 3 is the third-generation matter supermultiplet). The prediction 
of the top quark mass from this interaction was first considered by Ananthanarayan, 



Lazarides and Shafi fig . We find the picture which emerges from such an interaction 
to be very elegant. While the three Yukawa couplings A t) & iT are different at low energies, 
they evolve according to RG equations to a common unified value at the large mass 
scale of the grand unified theory. The scenario is reminiscent of the evolution of the 
three gauge coupling constants to a common value at the unification scale. While the 
gauge coupling unification leads to a prediction for the weak mixing angle, the Yukawa 
coupling unification leads to a prediction for the top quark mass. The top-bottom mass 
hierarchy then originates in the Higgs sector, so another prediction is a large ratio of 
Higgs VEVs. The weak mixing angle prediction has undergone several refinements as 
higher order corrections and a variety of threshold corrections have been considered. An 
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aim of the present paper is to compute such corrections to the top mass prediction to a 
similar level of accuracy. In particular we study: 

• the coupled two-loop RG equations for the three gauge couplings and the three 
Yukawa couplings. We give an analytic fit to the numerical results which is valid 
to better than 0.2%. 

• some implications of generating the top-bottom mass hierarchy through a large 
ratio of the VEVs of the two light Higgs doublets. This source for up-down mass 
hierarchy is generic in models which unify all three Yukawa couplings of the third 
generation. Since a priori there is no symmetry protecting the down-type Higgs 
VEV and consequently the down-type fermion masses, large radiative corrections 
to these masses typically arise and change the top mass prediction considerably. 
Such corrections will be suppressed if the squarks are much heavier than the hig- 
gsinos and gauginos; indeed, we identify two symmetries which could then protect 
the down-type VEV and masses. Whether such suppression is favored in models 
with large tan/3 is a question for future study [20|, and whether it is the case in 



nature will be determined by future experiments. 

two consequences of these large corrections to the b quark mass. For a certain range 
of MSSM parameters the m&/m T prediction cannot be brought into agreement with 
experiment; for a separate, smaller range, different GUT-scale boundary conditions 
must be used. 

the supersymmetric threshold corrections to the three gauge couplings and the 
three Yukawa couplings. These are given for an arbitrary spectrum of the su- 
perpartners of the minimal supersymmetric standard model, ignoring only the 
electroweak breaking effects in the spectrum. We find in particular that, when 
the above symmetries hold approximately, raising any or all of the superpartner 
masses increases rat for fixed a^{mz)- 

threshold corrections at the grand unified mass scale. We show that such correc- 
tions to the gauge couplings do not significantly affect the top mass prediction for 
a given a 3 (mz). We calculate the corrections to the Yukawa couplings from super- 
heavy splittings in the gauge 45, the 10 H and the 16 3 multiplets (these corrections 
are not very large), and give general expressions for further possible superheavy 
threshold corrections. 

the extent to which the predicted value of the top quark mass depends on assump- 
tion (III). This assumption is, in our opinion, the weakest part of the theoretical 
picture which underlies the top mass prediction (with the possible exception of 
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the electroweak symmetry breaking sector) . Even if the third generation only cou- 
ples through a single SO(10) invariant Yukawa interaction, the relation between 
this coupling and the values of A t , A& and A r renormalized at the unification scale 
may involve a set of mixing angles describing which two linear combinations of 
the doublets of the unified theory are light. An understanding of these mixing 
angles is in principle related to an understanding of why two doublets remain light 
(the doublet-triplet splitting problem). Our understanding of the resolution of this 
problem is at present not complete. We are however greatly encouraged by the 
following two facts: i) Due to the fixed point behavior of the RGE the prediction is 
somewhat insensitive to these mixing angles in a large class of models; ii) in very 
simple SO(10) models which provide a partial solution of the doublet-triplet split- 
ting problem, both Higgs doublets are in fact contained within the same irreducible 



representation (namely a 10 g ) |2T| , |20 



• the extraction of the b quark mass from experiment, using updated experimental 
information and including the dependence of the extracted value of this mass on the 
QCD coupling 0:3. This is of importance to us because the two crucial experimental 
inputs to the top mass prediction are nib and 0:3. 

In Sec. II we describe the basic framework for our calculation. A discussion of the 
implications of large tan f3 is given in Sec. Ill, where we examine two potentially very 
large corrections to the b quark mass prediction; however, if the MSSM parameters 
exhibit a certain hierarchy (for large tan/?) such corrections may be suppressed. In 
Sec. IV we give approximate analytic solutions to the two loop renormalization group 
equations, in the absence of threshold corrections. The extraction of the b quark mass 
from data is given in Sec. V. Armed with this experimental value of nib, we return in 
Sec. VI to the large corrections to the predicted nib, bounding these corrections (and 
consequently the MSSM parameters) and investigating the possibility of different GUT- 
scale initial conditions for the Yukawa couplings. The remaining threshold corrections 
are studied in Sec. VII, while in Sec. VIII the sensitivity of the predicted top quark mass 
to these threshold corrections is derived and discussed in a general way. In Sec. IX we 
give the prediction for the top quark pole mass and discuss its sensitivity to certain grand 
unified threshold corrections. In Sec. X we extend the discussion to include models in 
which the doublets are not completely contained in a single SO(10) irreducible multiplet, 
and to realistic models which include the other generations of matter. Sec. XI concludes. 
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II. FRAMEWORK 



The predicted value of the top mass depends on the renormalization group equations 
(RGE), the boundary conditions at the GUT and electroweak scales, and the threshold 
corrections at these two scales. The one-loop RGE in the minimal supersymmetric 
standard model (that is, with three generations, two Higgs doublets and no right-handed 
neutrinos) are given by 

16vr 2 ^ = £ K tv \l+ £ (la) 

0,1 v=t,b,T i=l,2,3 

i67r2 dlnAb = ^ K ^ l+ ^ W; (ib) 

at v=t,b,r i=l,2,3 

i67r2 rflnA I= ^ ^ + ^ ^ (ic) 

i/=t,b,T i=l,2,3 



and 



2 cHnj7i 2 



2 d\ng 2 



2 r/ln/ /:i , 



where £ = In p. and Jf t = (6, 1, 0),K b = (1, 6, 1), K T = (0, 3, 4), L t = (-§, -3 h> 



3 



L b = (— t^t, —3, — y), £ T = (— |, —3, 0), bi — ^ , 62 = 1, and 63 = —3. In our analysis, 
however, we employ two-loop evolution equations for both the Yukawa and the gauge 
couplings. The gauge couplings are related to experimental observables by g\ = |e 2 /(l — 
sin 2 ^), g\ = e 2 / sin 2 ^, and g\ = ^nas, whereas the Yukawa couplings are related 
to the running quark masses via A t = ^/2m t /vu, \b iT = \f2m b ^ T / and Vy + v 2 D = 
v 2 = (247 GeV) 2 . The VEVs of the two light Higgs doublets H v and H D are denoted, 
respectively, by Vjj and Vd, and their ratio is denoted by vu/vd = tan/? as usual. The 
boundary conditions in the gauge sector at the electroweak scale will be [p3jp4| the 
MS values sin 2 6 W = 0.2314 (appropriate in anticipation of a heavy top quark) and 
4:7i /e 2 = l/a em = 127.9; these are extracted from data using the 6-flavor standard 
model as the effective theory at the scale mz- In the Yukawa sector we have 



m T {m T ) = 1.777 GeV =>■ m T {mz) = 1.749 GeV, and the range of values derived 
below. The remaining uncertainties in sin 2 6w and a em , contribute negligibly to the 
prediction of m t . 

The tree-level initial conditions at the GUT scale Mq are g\ = g 2 = g% = gc f° r 
the gauge sector, but in the Yukawa sector they depend on the source of the low-energy 
Yukawa couplings. In SO(10) unification they typically arise from terms of the form 
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16 ? O 16 3 , where 16 3 is the chiral supermultiplet for the third generation and the Higgs 
multiplet O may be a 10 H or 126 H of SO(10). (The 120 H is antisymmetric and therefore 
makes no contribution.) To go beyond the SU(5) predictions and exploit the larger 
SO(10) symmetries, we must assume that both light Higgs doublets lie predominantly 
in a single SO(10) multiplet, rather than being arbitrary mixtures of doublets in several 
representations. As stated above, our understanding of which doublets remain light, 
and why they do so, is not complete. Thus to make progress we simply assume the 
mixing is negligible. Note that if the mixing is with other doublets which do not couple 
directly to the 16 3 , then the effect is only to split Xf from Xf T , and even a 30% splitting 
of this sort will have less than a 3% effect on the low-energy top mass, as we show 
below. Furthermore, a supersymmetric theory without mixing terms in the Lagrangian 
is technically natural, due to the nonrenormalization theorems. Then either O ~ 1Q H , 
in which case the initial conditions for the Yukawa couplings are 

At = Ab = A r = Ac , (3) 

or else O ~ 126 H from whence 3A t = 3Ab = A r = Xq. In the second numerical 
investigation shows that the prediction of m\ ) jm T is too low unless a^imz) > 0.13 and 
Trib(mb) < 4.0 GeV, as well as Ag > 4. This last requirement, however, when combined 
with the full SO(10) RGE, implies a Landau pole in the Yukawa coupling less than 20% 
above the unification scale, which in general allows 0{1) GUT-scale threshold correc- 
tions and makes any predictions impossible. In fact, the only option for perturbative 
unification with O ~ 126 H is for very large electroweak-scale threshold corrections to 
arise, which are of just the right magnitude and sign to restore the agreement of mfe/m T 
with experiment. Only for a very limited range of Ag values and MSSM parameters can 
such a scenario be successful, as we show in Sec. VI. For most of the parameter range, 
therefore, we need only consider the consequences of 16 3 10 H 16 3 . In this case, we must 
restrict Ag < 2 to make sure the Landau pole is at least a factor of 4 above the uni- 
fication mass; for higher values of Ag, we simply cannot make any reliable predictions, 
although these values are by no means ruled out. 

III. LARGE tan/3 

We discuss in this section some of the implications of the large value of tan/5 = 
vu/vd — 50 — 60 necessitated by the boundary condition Xf = X^ = X^. In the standard 
model and most of its extensions, the hierarchy m h <C m t (and similarly for the r mass) 
is imposed through Ab <C A t ; Ab is the small parameter quantifying the breakdown of the 
chiral symmetry which protects the bottom quark mass. In multi-Higgs models such as 
the MSSM, there is also the option of explaining m& <C m t by having the down-type Higgs 



7 



acquire a much smaller VEV than that of the up-type Higgs, vd <%■ We are forced to 
take this second route by our GUT boundary condition, which implies Ab ~ A t ~ 1. We 
will assume as usual that the electroweak symmetry is broken by an instability of the 
scalar potential m 2 v \E v ^ + m 2 D \H D \ 2 + Bfi(H v H D + h.c.) + |(# 2 + g' 2 )(\H u \ 2 - \H D \ 2 ) 2 in 
the Hjj direction. This breaking is then communicated to Hd by the soft SUSY-breaking 
term [iBHuH D : 

—aB , „ 1 Vn 1 

±sin2/3~ =_£~ (4) 



mfj + m 2 D 2 tan (3 vy 50 

So the magnitude of the hierarchy between the up and down sectors is determined by 
fiB /(mjj+m 2 -)), while its direction is set by the direction in which the instability develops 
in the scalar potential. The denominator is given by mfj + m 2 D = m\, the squared mass 
of the pseudoscalar neutral Higgs boson; evidently, either the \x or the B parameter or 
both must be much smaller than in order to generate the top-bottom mass hierarchy. 
How, and indeed whether, the various parameter values may arise will be discussed in a 



future paper [20 



At tree level, an appropriate choice of parameters in the scalar potential leads to 
vd <C vjj and hence to rrib,T m t . However, no symmetry has (yet) been imposed 
on the Lagrangian to protect such a hierarchy, and therefore we expect large radiative 
corrections. In fact corrections arise from the gluino- and higgsino-exchange diagrams 
of Fig. 1 (the analogue of Fig. la, where a bino is exchanged, is suppressed by the small 
hypercharge coupling). Such corrections have been overlooked in past work on large 
tan j3 scenarios, though they were discussed in the context of radiative mass generation 
[26|. Thus the MSSM prediction for becomes, after replacing Hjj — > vu/y/2, 

m b = A b ^ + 8mf + Smp (5) 



where 



'V2 



Sm (g) = ^nig Mb^= J r (m| >+ , mf _, m 2 ) (6) 



and 



5mi b h) = ^^ A * A t^| / ( m I+' m !-'/ i2 )' ( 7 ) 



I is given in the appendix, nig is the gluino mass, \i is the supersymmetric coupling of 
the two Higgs doublets, A t is the trilinear soft SUSY-breaking coupling of the stop fields 
to the up-type Higgs, and m^ ± and m t - ± are the squark mass eigenstates propagating 
in the loop. Numerically m g ~ Smu2 where mi/2 is the mass of either the gauginos at 
the GUT scale or of the wino at the electroweak scale. To appreciate the significance of 
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FIG. 1. The leading (finite) 1-loop MSSM contributions to the b quark mass. 

such corrections, consider the limit in which the squarks have roughly equal masses itlq, 
and /i or m g are either much less than or equal to m . Then I(ml, ml, 0) = l/ml and 
I (ml, ml, ml) = I /(2ml), anc ^ we nn d ; 



(Note that in the second line above, we have approximated tan/3 ~ 50, which is inac- 
curate if Smb becomes large and lowers the top mass prediction significantly; we will 
return to this point once we have extracted the experimental bounds on m&.) We see 
that the radiative corrections may in general be comparable to the tree-level mass; when 
we equate this prediction to the experimental value of m^ to extract Ab, we would find 
0(1) corrections and hence large changes in Ag and in the prediction for m t . Our final 
m t prediction would then be very sensitive to the exact values of the squark, higgsino 
and gaugino masses (and perhaps A t as well) — far more sensitive than expected from 
ordinary threshold corrections. On the other hand, the squarks may turn out to be 
relatively heavy, namely ml ^> [im- g , \iA t , in which case these corrections would be sup- 
pressed. For example, if m ~ 1 TeV but /i ~ m~ g ~ A t ~ 200 GeV then m\, changes by 
only ~ 6%, which in turn corrects the top mass prediction by ~ 4%. Of course, the sign 
of this correction is determined by the sign of the parameters which enter bm^. 

We should mention at this point that there is a diagram in the 2-Higgs standard model 
analogous to the higgsino-exchange diagram of Fig. 1, in which the stop propagator is 
replaced with a top propagator, the higgsino with the Higgs, and the couplings are 
replaced by A t \ t — > \ t and ji — > jj,B. For this diagram the large tan/5 enhancement 
gained by coupling the b to is manifestly and exactly cancelled by the fj,B/m 2 A factor 
from the propagators, independent of any symmetries. This threshold contribution is 
included in the function fn defined in Sec. V. 

We have seen that when ml ^> fim g , fiA t the large radiative corrections are sup- 
pressed. We also know that for large tan/3, m\ 3> fJ,B. These hierarchies can be 
related by imposing two approximate symmetries. The symmetry which sets n to zero 
is a Peccei-Quinn (VQ) symmetry under which the superfield H D and the SU(2)-singlet 
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bottom antiquark superfield b c have equal and opposite charges while all other fields are 
invariant. It is explicitly broken only by the soft SUSY-breaking term (iBHuHd and 
by the term nH v H D in the superpotential, so when treated as a spurion /i should be 
assigned a VQ charge opposite to that of H D . The SUSY-breaking term contributes 
at tree level to the VEV of Hp, while the supersymmetric one enters into (finite) loop 
diagrams which correct the b mass. We will quantify the degree to which this symmetry 
is broken by the dimensionless parameter e-pg = fi/mo. 

The symmetry which sets B to zero must also set the gaugino mass mi/ 2 and the 
trilinear scalar coupling A (and in particular A t ) to zero; note that mi/2 and A generate 
B through the RG evolution. The desired symmetry is in fact a continuous 1Z symmetry, 
and it is convenient to choose an 1Z such that the superpotential is invariant while the 
soft SUSY-breaking terms (except for the common scalar mass) is not. Furthermore, 
we choose an 1Z under which the scalar Hy is invariant but the scalar Hd and the b 
quark mass operator Qb c are not. We will assign the superspace coordinate 9 a charge 
of +1, the superfield Hjj a charge 0, and the quark doublet superfield Q a charge 0; 
then the superfield Hp carries a charge of +2, the superfield b c carries a charge of 0, 
and the superfield t c carries a charge of +2. Thus the left-handed quarks have charge 
— 1, the SU(2)-singlet b c antiquark has charge —1 and the SU(2)-singlet t c antiquark has 
charge +1. In a spurion analysis, the soft SUSY-breaking parameters m^, A, and B 
carry a charge of —2. Once again, we define a dimensionless symmetry-breaking measure 
t n = B/m < A/m ~ m 1/2 /m . 

When the VQ or 7Z symmetries are approximately valid, we find 

(f)(3)- (9) 

and 

, m f, „ ^ {eveeM (|) g) mt = A. (|) Q) mi . m 

As expected, the large tan f3 enhancement was cancelled by the e-pge^. factor, provided 
trig ~ A t ~ B and m\ ~ m\. 

If either of these symmetries is to hold even approximately, there must be a cer- 
tain hierarchy in the supersymmetric spectrum. Notice, however, that \x and 7711/2 are 
bounded below by LEP data |[27|| , at least for large tan/?: numerically, they read roughly 
(|/i| — m^/2) (|mi/2| — rnz/2) > niz/2, so we expect both \i and m.1/2 to be at least as 
large as rriz- Thus the VQ and 1Z symmetries can be meaningful only if the scalar 
superpartners are significantly more massive than the Z . This typically implies a degree 
of fine-tuning to get the proper Z mass, so it remains to be seen [[2(J whether it is more 
natural to expect m ^> rriz ~ fi ~ 1^1/2 ~ A t ~ B and fine-tune the Z mass, or to 
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expect mo ~ mz ~ /i ~ m^j ~ A t ^> B and fine-tune B. Perhaps electroweak sym- 
metry breaking does not arise from the usual running of the parameters in the Higgs 
sector, but rather from some other mechanism (which does not significantly alter the 
RG evolution of the gauge and Yukawa couplings). Ultimately, experiment will decide 
what if any hierarchy exists in these parameters. Some experimental information already 
exists: attaching a photon in all possible ways to the diagrams of Fig. 1 and inserting a 
flavor-changing vertex leads to an amplitude for b — > sj which again is proportional to 
e-pQe-ji tan f3. To reconcile this amplitude with the CLEO data on Y{b — > 57), we must 
either impose the VQ and 1Z symmetries or raise all the superpartner masses (since the 
operator for b — > 57 is of dimension higher than 4). We will leave these constraints 
to future work f20fl . (Analogous considerations can be extended to other observables 
such as B°B~°, D°Tf and K !^ mixing, the electric dipole moment of the neutron and 
proton decay.) For now, we can only make definite predictions of the top mass in the 
case TOq 3> [im- g , fiA t . 



IV. RUNNING AND MATCHING 

If we temporarily ignore all threshold corrections (the "unperturbed scenario"), the 
solution of the RG equations proceeds schematically as follows. By requiring that the 
two gauge couplings g\ and #2 meet, we solve Eqs. (|2"a| , |2T)| ) to obtain Mq and go- (At this 



stage there is only a weak dependence on the Yukawa couplings in the 2-loop RGE, and 
we may use representative values for them.) Then by running back down with Eq. ( ^c|) 
we predict Qzijnz) and hence a 3 (m z ). Next, we solve the two equations (|I^,|lcD for the 
two unknowns Aq and tan/3, using Mq and qq as well as m^intz) and m T (mz) and the 
initial condition (|3|). In practice, this step is simplified because tan/? will always be 
~ m t (mz) /mb(mz) ~ 50 — 60, so we may set sin/5 to unity from now on. We are then 
left with a single equation, namely m h {mz) /m T {mz) = R = Xbi^z)/ K{m>z), which we 
solve for the single unknown Ac- Finally, we use Ag to run down with Eq. (|Ta|) and 
obtain At(m^), which is then used to determine the top quark mass. 

More precisely, we will adopt the following procedure to correctly incorporate 2-loop 
RG evolution with 1-loop matching conditions. We choose to match the MSSM with the 
broken-electroweak standard model using as an intermediate step the 2-Higgs standard 
model (2HSM). This two-step procedure has two advantages: the presentation is clearer, 
and the various matching contributions are easy to isolate. 

• We first treat the gauge sector: we match the experimentally-determined gauge 
couplings gi{mz) of the standard model, which are essentially those in the 2HSM, 
to the couplings of the MSSM by integrating in the superpartners. The conversion 
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from the MS to the DR scheme is numerically insignificant. We then use 2-loop 
MSSM RGEs to run from mz to some GUT scale hg which we will fix to be some 
convenient value near 10 16 GeV; in this running we employ approximate values of 
the Yukawa couplings. We thus calculate the gauge coupling boundary values at 
this GUT scale. 

Starting with these gauge couplings and with a given set of Yukawa boundary 
conditions (collectively denoted by Ac for the moment) at the GUT scale [iq, 
we evolve the gauge and Yukawa couplings with 2-loop MSSM RGEs in the DR 
scheme to an arbitrary electroweak scale fxz and obtain Xt,b,r = Af fe s ^ M ' DR (/i^; Xq). 
These are then matched to the 2-Higgs standard model, in which the superpartners 
have been integrated out, to yield AjJ^. M ' DR (/i^; Xq, {n^a}) = ■^ S f M ' DR (/ i z; Aq)[1 + 
kt,b,r(nz] { m a})], where {m a } are the superpartner masses. 

We also start with the running MS value of the b quark mass m^ s (4.1 GeV), 
and evolve it with 2-loop QCD running [that is, in the no-Higgs low-energy stan- 
dard model (OHSM) having only strong and electromagnetic interactions] to ob- 
tain m° HSM ' MS (/zz) = m^ s (4.1 GeV) /rib- Similarly we run up the r mass to obtain 
^ohsm,ms = m T /r} T . Their ratio is defined to be R msM < MS (fi z ), which may 
be translated into the DR scheme: i? 01 ^' 55 ^) ^ R° nsM > m (nz) [1 - a 3 /37r]. To 
match the OHSM to the 2HSM requires some knowledge of the Higgs sector masses. 
At tree- level and in the limit of large tan f3, these are all given |28] in terms of mz, 



mw, and the mass mA of the physical neutral pseudoscalar Higgs: the (mostly 
up- type) scalar m^ ~ mz, the other (mostly down-type) scalar m# ~ mA, and 



the charged (also mostly down-type) scalars m H ± = y + m\. We then obtain 
R = R msM > m (vz; m A ) = R msM > m (^ z )[l + fniliz; m A )\- 

Next we demand R = Af SMM {fiz] X G , {m a })/ Xl RSum {fi z ; A G , {m a }) which we 
may solve for Aq. 

Finally, we use the A( HSM ' DR (/iz; Ac, {m a }) corresponding to this Aq and calculate 
the top pole mass m]? olc = [v/ v / 2)A^ HSM (/i^; A G , {m a }) [1 + ft{nz',rn A )), defined 
as the position of the pole in the 2HSM with perturbative QCD. The function f t 
contains the well-known contribution from perturbative QCD radiative corrections 
as well as often-neglected contributions from Yukawa radiative corrections. Note 
that the observable m]? ole must be independent of \iz to 1-loop order; we have 
indeed checked that our final values do not change by more than a GeV or so 
when we vary \iz around the electroweak scale. To be specific, we will use the 
value Hz = m z in all the explicit values we present below. 
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We will make the following approximations when appropriate. First, we use the full 
1-loop threshold expressions involving A^, A^, A^, and grf, with the following exception: 
when integrating out the superpartners and matching to the 2HSM, we are neglecting 
operators of dimension > 4 which are suppressed by the superpartner masses; when 
calculating the top (rather than bottom or tau) mass, this amounts to neglecting finite 
terms of order n|/ , m.superpartner> which is not valid in the case of a light superpartner but 
will numerically be sufficiently accurate. Second, corrections proportional to electroweak 
gauge couplings have only been included in leading log(msusY) approximation; this 
means that we have neglected finite parts from SUSY loops and all electroweak gauge- 
boson loops. Third, we keep only the dominant diagrams (namely the ones in Fig. 1) 
from the class of finite diagrams proportional to e-pQe-%, €p Q and e\\ the rest contribute 
at most ~ 1% effects, or much less if e-pq, < 1 and en < 1. Fourth, when integrating 
out the Higgs sector we neglect the effects of various 1/ tan/3 mixings. Finally, in the 
numerical results we have grouped together and assigned common average masses to the 
squarks (rriq), the sleptons (m|), and the higgsinos (m^). 

The evolution of the Yukawa and gauge couplings from \i G to fiz in the MSSM, even 
to 1-loop accuracy, cannot in general be calculated analytically. Numerical solutions are 
straightforward, and show that, in order to obtain the correct m b /m T ratio, the GUT- 
scale Yukawa couplings must be 0(1), so the top mass is typically predicted near its fixed 
point value fl6|]T7l]24jl . With this in mind, it will be useful and illuminating to obtain 
simple, explicit approximate solutions by fitting Xt^A^z) numerically to a quadratic 
polynomial in 1/Aq. The constant term reflects the independence of the low-energy 
Yukawa couplings on the initial conditions in the large- Aq limit, a consequence of the 
fixed point behavior. When varying Aq between 0.5 and 2, we find in the unperturbed 
scenario 

XtfiAvz) - A,b, T + ^v^ 1 + (11) 

where A t>b>r = (1.099,1.014,0.673), B t)bjT = (-0.045,-0.012,-0.107), and C t>b>T = 
(—0.019,-0.025,0.001); these values correspond to fiz = 90 GeV and to g% = go 
(that is, no superheavy threshold corrections) which leads to a 3 (m^) = 0.125. We 
have checked that the errors we make in this fit are smaller than 0.2% over the 
entire range 0.5 < Ag < 2.0. We may then solve the quadratic equation R = 
(A b + B b /X G + C b /\Q)/(A T + B T /\ G + C t /\q) to obtain Ac as a function of R, which 
in turn gives an explicit solution for A t (yUz) as a function F(R) of the experimentally- 
determined ratio R. The functions (|i~l~D, as well as the ratio R = \{^z) I 'K(nz), are 
plotted versus Ac in Fig. 2a, while F(R) is shown in Fig. 3a. 
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FIG. 3. The induced dependence \t( m z) = F(R) derived from Fig. 2, as well as 
the sensitivity functions {Fn\R)} and {F„ k \R)}. 
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V. THE MASS OF THE BOTTOM QUARK 



What is the experimental value of BP. Using updated values from the Particle Data 
Group compilation |[29|| , we have reanalyzed the well-known extraction of the b mass 
from e + e~ collisions via QCD sum rules p0| . The idea is to use a dispersion relation to 
relate the experimental spectral distribution of e + e~ — > bb to the expectation value of the 
product of two vector 6-quark currents. This product is then rewritten as an operator 
product expansion; perturbative QCD is used to calculate the coefficients of the identity 
and other operators in the expansion, while the nonperturbative information is assumed 
to be contained in the condensates of these other operators. For the b system, one 
expects the nonperturbative terms to be at most 0((ipip) /ml ~ Aq CD /??7|) and therefore 
negligible. The remaining calculation is purely perturbative QCD, so the uncertainty 
in nib will be dominated by our ignorance of the 0(a 2 ) terms in the calculation of the 
coefficient of the identity operator. This coefficient can be calculated reliably only for 
highly off-shell momenta, for instance q 2 <C m 2 . Expanding the coefficient in powers 
of q 2 /m 2 results in a relation between the moments of the spectral distribution and 
derivatives of the coefficient at q 2 = 0: 
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where the sum approximating the spectral integral is over the various resonances V 
characterized by a mass My and a width IV, Et — 10.56 GeV is the estimated threshold 
energy where the continuum bb production begins, and M.o n and A n are numerical 
constants given in reference |50| . 

The mass is extracted from the ratios of the first few successive moments, 
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(14) 



where r 0j2 = -0.00452, r 0j3 = -0.00462, Ql = -0.0286, a 2 = -0.197, and k has not yet 
been calculated but is expected to be at most 0(1) for the first few moments. We will 
make a very rough estimate of our uncertainty by allowing k to vary between —2 and 2. 
The strong coupling a 3 in these expressions must be run down from its given value at 



mz- The parameter is the Euclidean pole mass of Georgi and Politzer [[H], which 
is related to the running MS mass we need via 



m btE = m b (fi) 
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Using the first three moments, we obtain 3.93 GeV < m&(4.1 GeV) < 4.36 GeV if 
a 3 (m z ) = 0.11, and 3.86 GeV < m b (4.1 GeV) < 4.42 GeV if a 3 (m z ) = 0.12. The 
central values of m{,(4.1GeV) extracted from the next few moments fall well within 
these ranges, providing some confidence that our error bars are not too small. Thus we 
estimate 



m b (m b ) 



4.15 GeV ±0.22 GeV, a 3 (m z ) = 0.11 
4.14 GeV ±0.28 GeV, a 3 (m z ) = 0.12 



(16) 



The central values we extract for m;, are not very sensitive to a 3 since the a n a 3 <C 1; 
they are somewhat lower than in the older analyses mainly because the more precise 
experimental value we use for the electronic partial width of the T(9460) is higher than 



the older value. Our error bars in m& are larger than those of Gasser and Leutwyler [32 
because of the different ways we estimate the error from the 0(a\) terms. 

The 2-loop QCD evolution between 4.1 GeV and \x z reduces the mass of the b by 
a factor r\ h ~ 1.437 + 0.075[a 3 (m z ) — 0.115]/0.01 for \i z = 90 GeV; the corresponding 
electromagnetic reduction factor for m T is r] T ~ 1.016. The translation from MS to 



DR increases rrib by roughly half a percent and has virtually no effect on m T . Finally, 
we match the 0HSM model to the 2HSM by including the radiative corrections of the 
Yukawa couplings via the function fR(m A ), 



/i?(mA ) =jR 2HSM / ^0HSM_ 1 

~ — 0.014 \n(m A /m z ) when \x z = m z and m A > va z 



(17) 



and arrive at the final value of R: 
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' 1.67 ±0.09, 


a 3 (m z ) 


= 0.11 


1.58 ±0.11, 


a 3 (m z ) 
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1.62 ±0.09, 


a 3 (m z ) 


= 0.11 


1.53 ±0.11, 


a 3 {m z ) 


= 0.12 



[1 + f R {m A )\ 
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1000 GeV 
1000 GeV. 



The exact expression for may be found in the appendix. Notice that a heavy second 
Higgs decreases the apparent experimental value of R, or more intuitively increases the 
GUT prediction of R — hence to agree with the experimental value, we need to lower the 
prediction of R, which entails raising X G and with it m t . The lightest top masses result 
from a light Higgs sector. 
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VI. IMPLICATIONS OF THE LARGE CORRECTIONS 



With the experimental value of rrib in hand, we can return to the threshold corrections 
of Eq. (||) to bound the allowed range of MSSM parameters and to possibly allow different 
initial conditions for the Yukawa couplings at the GUT scale. 

First, let us see how large can the corrections to mj become before m;, as pre- 
dicted from Yukawa unification disagrees with the experimental range given above. 
We focus on the (usually dominant) gluino contribution; similar considerations will 
bound the higgsino diagram. It is convenient to first divide Eq. (Q) by m T , and sub- 
stitute tan/5 = (m t /m T )(X T /Xt) = (m t /m T )(X h / X t )RussM - -^Russu( m t/ m r) and 
H m l + 5 m f _) m §) = l/ m cfr> to obtain: 

p m b p , f n nr ,8 a 3 \ m t (R MSS M) (m s n\ 

-Kexp — — — -Kmssm + u.yo-— — 5— 

m T V SAirJ m T \ m is J 



m t (R MS sM) ( mgfi\ 

-Kmssm H 7 „ p — — — 

73GeV V m cfr/ 



(19) 



Note that m t depends on X G which is in turn determined by -Rmssm- The R exp on the 
left-hand side is the experimental value R extracted above, while -Rmssm on the right- 
hand side denotes the value of Ab/A T obtained by running down in the MSSM from the 
GUT scale. Also, to first approximation m t (i?MSSM) — (174 GeV) At (.Rmssm)- As we 
will establish from Eq. ([29 bp below, and as illustrated in Fig. 2a, -Rmssm is bounded 
from below by roughly 1.6, corresponding to Ag — > 00; and from above by roughly 
2.4, corresponding to Ag — > 0. Since the latter limit also corresponds to m t — > 0, it 
can be improved by enforcing the experimental bound m t > 130 GeV, so we conclude 
that 1.6 < -Rmssm < 2.15. We use these bounds and those of Eq. (|I8D to set limits on 
(rngfi/m^g): if this quantity is positive, an upper bound results from taking the smallest 
-Rmssm possible and the largest -R exp allowed, in which case m t is fixed at its maximal 
value; if that quantity is negative, the largest .Rmssm and smallest -R exp are needed, in 
which case mt — 130 GeV. In this way we find the surprisingly stringent limits 

-0.37<^< 0.08 (20) 



m 



cff 



which are obviously phenomenologically interesting signatures, but are also relevant to 



the electroweak symmetry-breaking sector of this model [[20| . 

Second, recall that if we impose the GUT-scale initial condition 3A t = 3Ab = A r = Ag 
corresponding to mass generation through a 126j j, then after evolving down to the 
electroweak scale the prediction for R is too small. Large Srrib corrections can restore 
agreement with experiment. Fig. 2b shows .R as well as A t , Ab and A T at the electroweak 
scale as functions of Ag- In the presence of the 126 H , the coupling Ag must be kept below 
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~ 1.5 (rather than ~ 2 for the 10 H ) to raise the Landau pole by a factor of 4 above 
the unification mass. Also, since m t ~ (174GeV)A t , the lower bound of m t > 130 GeV 
implies A t > 0.75 and hence from the figure Xq > 0.82. Within this restricted range, 
-Rmssm varies between 1.0 and 1.1, so to reconcile this with i? exp (now using Ab/A t — 0.89 
and again focusing on the gluino diagram) requires 

0.24<^<0.42. (21) 

™cff 

We learn that if the gluino, higgsino and squark mass parameters satisfy mj/j/m^ = 
0.33 ± 0.09 and if A G — 1.1 ± 0.3 then the mass may originate perturbatively from 
the coupling 16 3 126 H 16 3 . Notice that in such a scenario, even if mjju/m^ is known 
precisely, then the uncertainty in R exp is large enough that the top mass prediction will 
usually be very imprecise. 

Thus for Yukawa unification in the MSSM it is useful to distinguish four regions of 
parameter space: 

1. If \rrigiJ, I 'ml s \ <C 1 then the Srrib corrections may be ignored, the mass must orig- 
inate from a 16 3 10 H 16 3 interaction, and we can predict m t with little further 
dependence on the MSSM parameters, as shown below; 

2. If —0.37 < mgfx/ml s < 0.08 then the mass must still arise from a 10 H , but now 
the prediction for m t depends very sensitively upon mgfi/m^Q and can vary over 
the full experimentally-allowed range; 

3. If 0.24 < m~ g /i/m 2 cS < 0.42 then the 126 H must be used, with Ag — 1.1 ± 0.3. 
In this case the prediction, while imprecise, tends to be in the lower half of the 
experimentally-allowed range; 

4. If rrigfi/mlff lies outside of these ranges then perturbative Yukawa unification under 
our assumptions cannot be reconciled with experiment. 

VII. THRESHOLD CORRECTIONS 

Next, we investigate the deviations in the top mass prediction induced by threshold 
corrections at the GUT and SUSY scales. For convenience, we choose to always match 
the full SO(10) theory with the MSSM at the scale Mq of the unperturbed scenario, 
where g\ and Q2 met: we define from now on fie = 2.7 x 10 16 GeV. The top mass 
prediction as a function of R is completely determined once we know the initial conditions 
{A^, A^, X^, gf, g^i df} a t the fixed scale fie, as well as the functions ^^({ma}) °f the 
superpartner masses which match between the MSSM (in which we run with the RGE) 
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and the 2-Higgs standard model (in which we calculate the top pole mass). So we first 
calculate the perturbations et t b,r,i,2,3 to the initial conditions and the matching functions 
k tybyT in terms of the various mass thresholds and a^{mz)- Then we study the sensitivity 
of the top mass prediction to these e's and fc's. A linear analysis of these perturbations 
is sufficient for our purposes. 
Define for convenience 

t a = \a(m a /m z ), T a = ln(M a /// G ), (22a) 
t a ,b = G x {m 2 J n\,m\l n\) 

ln(m a /// z ) - |, m a > m b , 

\n(m a=b /u z ), m a = m b , (22b) 
]n(m b /fiz) ~ f , rn a < m 6 , 

*' oAc = In 1 " »' ^ (22c) 

Hz 

where G\ is given in the appendix, {m a } are the masses of the various superpartners 
and {M a } are those of the superheavy particles. The functions t a and t a ^ b yield the 
exact threshold corrections and will be used in the dominant terms in k t b T below. The 
function t' a b c only yields threshold corrections in the leading log(m SU sY) approximation 
and will be used in the subdominant terms. The contributions of the superpartners to 
the gauge /3-functions will be denoted by b\. Superheavy threshold corrections can arise 
from couplings dressing the 16 3 10 g 16 3 vertex and having strength gc (the gauge 45). 
Xg (the 16 3 and 10 H ), or some other \' A . We write the corresponding contributions 
of any superheavy particle with mass M a to the Yukawa RGE by L^g G , K^Xq and 
HaK'v A \% We denote squarks, sleptons, higgsinos, winos, binos and gluinos by q, £, 
h, w, b, and g, and define olq = QgI^ anc i Vx = ^l/^ f° r an y x - We will neglect all 
electroweak-breaking effects in the mass splittings, and the few subdominant corrections 
mentioned below. We expect all such effects to alter wit by less than a GeV or so. 

The gauge couplings at fie may be completely determined by their low-energy values 
0-i,2,2,{ m z) and by the SUSY spectrum simply by running them up from the Z mass to 
fie- When the superpartner masses are changed away from the 90 GeV value of the 
unperturbed scenario, the couplings gf and gf at fie deviate from the common value 
gc = 0.730 of the unperturbed scenario in an easily-calculable way. Similarly, gf deviates 
from this go value when the masses of the colored superpartners changes and also when 
we vary |33| a 3 (mz) away from 0.125. While these changes to gf are calculated using 



the experimentally-accessible quantities {m a } and a 3 (m^), from a top-down viewpoint 
they should be regarded as the net threshold corrections resulting from integrating out 
the superheavy (or Planck-scale []34]]) degrees of freedom in the SO(10) theory to arrive 
at the MSSM. The initial conditions for the gauge couplings become 
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9iG = 9g{^ + ei) = 9g 
92G = flte(l + £2) = 9G 
93G = 9g(1+£3) = 9g 



, a G /ll 9 2 
1 - -r 1 t^s + + -=t 
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47T 
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2« 3 K) V "3 



(23a) 
(23b) 
(23c) 



where 5^3/03 = [0(3 (m#) — 0.125] /0.125. (Note that the fermionic superpartners w, h 
and g must be integrated out at \/2 times their mass.) The Yukawa couplings at Mq 
differ from Ac due to GUT thresholds effects, and hence 



Kg = Ag(1 + e t ) = A G 
Kg = Ag(1 + £5) = Ag 
Kg = Ag(1 + e r ) = A G 



1 + i E Wvo + L " ao + E 

1 + TZ + ^T a G + E 



(24a) 
(24b) 
(24c) 



Two remarks should be made at this point. First, the SUSY threshold corrections 
actually enter Eqs. ([24]) indirectly, even if all superheavy particles are degenerate, since 
these corrections generically change the scale at which g\ and g 2 meet and therefore 
change the predicted average mass of the superheavy particles. The result is a shift in 
all the T a that is independent of a. Second, any effect that contributes equally to e t , e& 
and e T makes no contribution to our final result, since this only amounts to a redefinition 
of Ag- 

Finally, the superpartner masses induce threshold corrections to the Yukawa cou- 
plings when they are integrated out of the MSSM to yield the 2-Higgs standard model. 
We find 

4nk t = y t [\t~ tR ~ h + t hJl ) + \y b t- bR ; h + |a 3 (t h ~ g + £ t ~^) + 3a 2 fa^ + ft^- - f^) 

(25a) 



_U 3 n, 1+ - 4- i6 + -4-9+ 16+' i V 
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4vrA; T = y r (§t fj? ^ + ^ ^ + 3a 2 (fa,* + \t K , A - t'- -^ 



(25b) 



(25c) 

where the y v and a% are evaluated at /i^-The non-logarithmic threshold correction of 
Eq. (H) must be included if the squarks are not much heavier than /i, trig and A t : 



tan /3 ( 8 m g /j A 2 
47r I 3 3 mL 47r m 2 



(26) 
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where = l//(m| , rn~_, m|) and m e 2 ff , = l//(m~ + , m|_, /x 2 ). 

Note that the majority of the terms in k v are due to wavefunction renormalization 
from scalar-fermion interactions and hence increase with the superpartner masses; the 
eventual conclusion will be that the smallest m t is predicted when all superpartners are 
as light as possible (if k' b may be neglected, as discussed in Sec. III). 

The mass of the Higgs bosons (which are determined by mw and tua) enter the top 
prediction through the matching of R between the OHSM and the 2HSM, and through 
the calculation of the top quark pole mass in the 2HSM. The former was included in the 
previous section as a correction to R, while the latter is studied in Sec. IX. 



VIII. SENSITIVITY FUNCTIONS 



When the initial conditions are perturbed, so are the coefficients of the fit in Eq. ([TTJ) . 
We expect the fit parameters to vary linearly with the perturbations as long as these 
are sufficiently small, and so we write 



(A v + ^Ojy n e n) H 7 1" , > 



;i + k) (27) 



where the sum ranges over n = t,b, r, 1, 2, 3. The new fit coefficients {a un , b un , c un } must 
be computed numerically. We then solve R = \ h (m t ) / X T (m t ) for Ag as before, substitute 
back into Eq. ( p?[ ) and expand to first order in e n and k v to find 

At(/xz) =F(R)+J2Fi e) (R)e n + J2^ k) (R)k v . (28) 

n v 

F(R) and the five "sensitivity functions" F^ Tl 2 (R) are plotted in Fig. 3a while the four 
sensitivity functions F^ e \R) and F^ T (R) appear in Fig. 3b. We have checked that, for 
the entire range 0.105 < a^(mz) < 0.13 and {m a } < 3 TeV, our approximation is off by 
at most ~ 1%. 

We can understand the behavior of A t (/x^) shown in Figs. 2a and 3a,b as follows. 
• First, as is well known ||15|| , A t and Ab (and to a lesser extent A T ) are quite insensitive 



to Ag for large Ag, since they both tend towards a fixed point as Ag — > oo. The 
fixed-point behavior of A t is manifested in the smallness of F} e \R) — changing 
the initial Xf by 10% changes its final value A t (/iz) by at most 1%. The sensitivity 
of A t (/Uz) to A^ is even less for small R, since small R implies large Ag and hence 
a stronger fixed-point behavior for A t . 

The value at the fixed point is determined predominantly by at low energies, 
and hence F% (R) is large; conversely, aip are smaller at low energies and hence 
have little influence, so F[ e l(R) ^ 1. 
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Next, we establish the fixed-point nature of Xb(t)/X T (t) = R(t) by examining its 
RGE: 

i6 ^ 2 rfln_R(t) = A , + _ + _ 16 g , (2Qa) 

~ A? ( 4 - JO - f al (29b) 

where in the last line we have made the rough approximations A t — Ab and g\ ^> g\. 
For very small A G the first term is negligible, and R(t) is driven purely by the gauge 
coupling evolution from R G = 1 to R(/iz) — 2.4, independent of Ac- For very 
large A G , the first term dominates at the beginning, but now R decreases almost 



instantly until R ~ y 3/4 while A t decreases quickly and becomes independent of 
A G ; the subsequent evolution from these effective initial conditions to R(^z) — 1-6 
is therefore also independent of A G . Thus R(t) has a fixed-point behavior as a 
function of A G for both small and large A G . The values of interest to us, 0.5 < 
A G < 2.0, lie in the intermediate- to large- A G range, which is why R(nz) becomes 
less sensitive to A G as A G — > 2. 

• To understand F(R), observe that both R and A t display similar fixed-point de- 
pendence on A G , so when we eliminate A G to obtain X t (/iz) = F(R) we arrive at a 
roughly linear dependence of A t on R. 

• The behavior of F^J(R) follows from similar reasoning. For small A G , changing R G 
via efc — e T has a reasonably large effect on R(m t ), which requires (in order to match 
the experimental value) a moderate change in A G — see Fig. 2a. For large A G , the 
same change in R G has only a small effect on R(nz) due to the fixed-point behavior, 
but now this small effect gets magnified back into a moderate change required in 
A G . Thus a fixed — e T always necessitates roughly the same change in A G . But 
the resulting change in X t (nz) is tiny for large A G or equivalently for small R, and 
that is why F$(R) become small for small R. Now F t (e) (R) + F^ (R) + (R) = 
because when t t = q, = e T = e the final value of mt must not depend on e (this just 
amounts to a redefinition of A G ). Hence (R) = -F^{R)-F t (e) {R) ~ —F^ e \R). 

• Finally, F ( t k) is by definition equal to F(R), while F ( h k \r) and F$ k \R) measure the 
changes in A t (/xz) induced by changing R directly at the top mass scale, so their 
behavior follows immediately from the dependence of A t on R described above. 

IX. PREDICTIONS FOR THE TOP MASS 

The final step is the calculation of the position of the pole in the propagator of the top 
quark within the 2-Higgs standard model. (Note that the pole mass is only a parameter 
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in the calculation of experimental observables; we leave the study of the relation between 
m^ ole and the actual "top mass" extracted from collider data to future work.) There are 
two important radiative corrections to the pole mass: the usual QCD correction from 
gluon dressing, and Yukawa interaction corrections to the top quark propagator and to 
the Fermi constant. The result may be written as 

+ ^ + R + (1 + ^QCd) 

= A t 177GeV [1 + f t (m A )} . (30) 

In the last line we have substituted A^ HSM = A t and 1 + 5qcd = 1 + 5a3/3n — 
(8a 3 /47r) \n(m t /fi z ) — 1.015 when \i z = m z , and defined /t(m^) = T lW (0)/2m 2 v + 
8m t /m t + \5 l L + One should not confuse the pole mass in the above equation 

and the Euclidean pole mass defined in Eq. (15): the one just above represents the 
real pole of the propagator at timelike momenta calculated in perturbative QCD, while 
the Euclidean pole mass does not actually correspond to any pole in the propagator 
(for obvious reasons). Notice also that the above top pole mass, which is scheme in- 
dependent, has been written in terms of DR quantities. The dominant 1-loop cor- 
rection to the muon decay constant (used to define Gp) is taken into account via 
247 GeV = l/foWGp = v [1 - E w (0)/2m^] where E w {0) is the top contribution 
to the self energy of the W at zero momentum. The wavefunction renormalization 
of the top quark propagator is given by and while the mass renormaliza- 
tion 8m t is /i^-independent and vanishes in the 't Hooft-Feynman gauge we employ. 
The complete expression for f t is given in the appendix; to a good approximation, 
ft ~ 0.04 + 0.003 In m A /fi z . 

The complete form of our prediction for mt may be obtained by substituting the 
expressions for R, e n , and k„ into Eq. (p8|) and inserting the result into Eq. (|30"|). To 
untangle the various dependences of the prediction on the SUSY- and GUT-scale param- 
eters, we divide the various contributions into three classes, and discuss each separately 
before combining them into a prediction. 

First, let us ignore both the (potentially large) finite corrections discussed in 
Sees. Ill and VI as well as GUT-scale thresholds, and concentrate on the logarithmic 
SUSY-scale threshold corrections. Varying only the superpartner masses affects the 
predictions in three ways: through k„, through e$, and (as remarked above) through the 
a-independent shift in all the T a . The latter effect (which actually shifts all GUT-scale 
masses together) is shown below to be small. The first two can be significant, but only 
when they serve to increase m t . For example, if a^,{m z ) = 0.11, mb(mb) = 4.2 GeV 
and m>A = 90 GeV, then when all superpartners have a mass equal to m z we predict 
m t = 167 GeV. By allowing the various {m a } to vary between m z and 3 TeV, we obtain 
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m t as high as 195 GeV but only as low as 163 GeV. To understand this fact, recall that 



F^ e \R) ^> F^'iR) ~ Fi £> (R) so we may neglect relative to e 3 . In the three k v the 
dominant terms are those proportional to y t , Ub and 0:3. Therefore, for fixed a 3 (mz) the 
top mass prediction depends predominantly (to within a few GeV) on the squark, gluino 
and higgsino masses rather than the slepton, wino and bino masses. Keeping only these 
terms, we find the approximate formula 

5a 3 



K — F + 



2 4vr 



l F ( e )«G 

2 3 a 3 



IF. 



(e)«G 
47T 



+ t, 



-^(F + Fi kh 
34vr 



(3F + F b (k> 



2 4tt V 



3F b (k) ) 



(31) 



(k) 

where the R dependence is implicit and we have used = F. As it turns out, the 
positive terms proportional to t^g and t-^ (from wavefunction renormalization in the 
k u SUSY threshold corrections) are always larger in magnitude than the negative terms 
proportional to tg and t^ g (from e 3 SUSY threshold corrections to the QCD gauge 
coupling). Thus — if the 6m b corrections are ignored — heavy superpartners inevitably 
lead to a heavier top. 

by a 3 (mz), which changes not only 



The top mass prediction is further influenced 
the value of g 3 used in the running of the Yukawa couplings but also the value of R as a 
function of m b . Finally, the Higgs mass parameter enters into the matching between 
R in the low-energy no-Higgs standard model and the R in the 2-Higgs standard model, 
and into the expression for the pole mass of the top quark. The resulting dependence of 
m\ ole on nig, m g , /i ~ m^, m^, m b and a 3 (m z ) is somewhat complicated, so we display 
it in two complementary ways. First, we list in Table I the predictions of our complete 
expressions for various choices of the parameters, again omitting the 5m b corrections of 
Sec. III. (Such an omission is clearly unjustified for many of the parameters chosen for 
this table — it is only intended to illustrate the conclusions of the previous paragraph.) 



In this table we have used as average values 



mi 



mi 



500 GeV. Second, we 



(k) iii 

can approximate F, F b and F 3 as linear functions of R, and further approximate 
y t ~ F 2 /4ti and y b ~ 0.87y t . For a given set of values for {m^, a 3 (mz), m b } we obtain 
an expression 
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m t 
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c 3 (t q 
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V2g, 
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+ C q,h^q,h- 



(32) 



The values of the q are given in Table II, and the resulting top mass values are accurate 
to within ~ ±5 GeV when the various masses are varied between mz and 3 TeV. 

Next, we consider the finite corrections to the b quark mass as discussed in Sec. III. As 
we saw, these are small only when the squared squark masses are much greater than m g fi 
and \iA t . In this case, choosing for concreteness m~ h ~ \i = 100 GeV, m g = 300 GeV, 



100 GeV, mo- 



rn; 



1000 GeV and = 1000 GeV, and considering several 
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TABLE I. The top quark pole mass predictions for a range of values of the parameters. 
All masses are in GeV, and we have used set m,^ = m~ h = nig = 500 GeV, though the results 
are almost independent of these masses. We have not included the 5m b corrections which are 
large in some cases considered in this table. 



m A 
a 3 (m z ) 
m b {m b ) 




90 GeV 

0.11 
4.4 GeV 


1 TeV 
0.11 

4.4 GeV 


90 GeV 

0.12 
4.4 GeV 


1 TeV 
0.12 
4.4 GeV 


90 GeV 

0.11 
4.0 GeV 


1 TeV 
0.11 
4.0 GeV 


90 GeV 

0.12 
4.2 GeV 


1 TeV 
0.12 
4.2 GeV 


nig nig 


m h 


















90 90 


90 


157 


164 


174 


179 


172 


177 


180 


184 


90 90 


3000 


166 


172 


183 


188 


181 


185 


188 


( a ) 


90 3000 


90 


161 


168 


181 


186 


177 


182 


186 


( a ) 


90 3000 


3000 


169 


177 


189 


195 


185 


190 


( a ) 


( a ) 


3000 90 


90 


173 


181 


194 


( a ) 


189 


( a ) 


( a ) 


( a ) 


3000 90 


3000 


177 


184 


197 


( a ) 


193 


( a ) 


( a ) 


( a ) 


3000 3000 


90 


158 


167 


181 


187 


178 


183 


188 


193 


3000 3000 


3000 


162 


171 


184 


190 


181 


186 


191 


( a ) 



a For these parameter values, Ag > 2. 



TABLE II. Coefficients of the approximate formula for the top mass, Eq. (]32]), which is 



accurate to ~ ±5 GeV when the various masses are varied between niz and 3 TeV. 



m A 

a 3 (mz) 
m b (m b ) 


90 GeV 

0.11 
4.4 GeV 


1 TeV 
0.11 
4.4 GeV 


90 GeV 

0.12 
4.4 GeV 


1 TeV 
0.12 
4.4 GeV 


90 GeV 

0.11 
4.0 GeV 


1 TeV 
0.11 
4.0 GeV 


90 GeV 

0.12 
4.2 a GeV 


1 TeV 
0.12 

4.2 a GeV 


m? ( GeV) 
c 3 (GeV) 

^(GeV) 


154 
4.8 
8.6 
3.9 


162 
4.3 
8.0 
3.8 


173 
4.2 
7.9 
3.8 


179 
3.8 
7.4 
3.8 


170 
3.8 
7.4 
3.8 


177 
3.4 
6.9 
3.7 


179 
3.8 
7.4 
3.8 


185 
3.3 
6.9 
3.7 



a For m b (m b ) = 4.0 GeV and ot 3 (mz) = 0.12, most parameter choices lead to Ag > 2. 
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values of 0:3 (toz), we obtain the predictions shown in Fig. 4 as solid lines. The pole 
mass of the top quark is shown as a function of the running mj parameter (discussed 
in Sec. V) indicated on the upper horizontal axis. The allowed range for mi, according 
to Eq. (|I8D , using a 3 (mx) = 0.12, is also shown on this axis. We learn that if the 
SUSY parameters are sufficiently hierarchical that 5mf, corrections may be neglected, 
then the top quark is predicted to be heavier than ~ 175 GeV for a^(mz) > 0.115. 
It should be remembered that the prediction for a^{mz) without GUT-scale threshold 
corrections is around 0.124, so values much lower than this correspond to large GUT-scale 
corrections to g%. Furthermore, we find that perturbative Yukawa unification demands 
«3 < 0.125. These last two observations are in rough agreement with previous authors 
24] , |35[| . If it turns out that 5m& is significant, our predictions may change considerably 



and become highly dependent on the SUSY parameters. If 5m b > 0, then k' b > 0, and 
since F b k \R) > the top mass prediction can only increase; that is, either the change 
is small and the top stays near its maximal value of ~ 180 — 190 GeV, or else the change 
is too big and the corresponding SUSY parameters are excluded (see Sec. VI). On the 
other hand, if Srni, < the top mass prediction can be significantly reduced. We show 
in Fig. 4 the predictions that result from a light squark spectrum, namely m~ h ~ /i = 
250 GeV, m-g = 300 GeV, = 100 GeV, m~ q = m i = 400 GeV and m A = 400 GeV. 
The appropriate horizontal axis is now the lower one, which is obtained from the upper 
axis by — > m b + bm b . (Since 5uib depends on tan/3 and thereby on m t , we use 
the central prediction of m t in the figure as a rough guide in rescaling the horizontal 
axis. Also, we have checked that the fit in Eq. ( pT|) is still reasonably valid for these 
low values of m t .) As is evident from the rescaled bounds on m;, shown in the figure, 
the top mass prediction now would encompass essentially all the experimentally allowed 
range. In other words, we obtain a meaningful prediction only if the squarks are much 
heavier than m^/i and [iA t or if we know that bmb > 0. We will discuss whether these 
requiremnts are likely to be satisfied elsewhere |2^ . 

It is important to note that our prediction for m t , for given experimental inputs and 
for Snib ~ 0, is considerably larger (by 10 to 20 GeV) than has been previously obtained 
using lowest order analyses ||TB|,|TTJ. Carena, Pokorski and Wagner [^j have briefly 
considered the condition A t = Ab = A r as a particular case of GUT boundary conditions, 
and employed 2-loop RGEs to numerically obtain top mass values with which we agree 
in the minimal scenario. They do not, however, describe the complete dependences on 
nib and a% nor do they address the question of 5mb or of superheavy corrections. Finally, 
we reach different conclusions about the dependence on the superpartner spectrum. 

We turn now to the threshold corrections which may be present at the GUT scale. 
These fall into three classes, corresponding to the three terms that make up e t ^, T in 
Eqs. (0): the splitting of the 16 3 and of the 10^ which contribute in proportion to 
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FIG. 4. Our predictions for mP t ° c without superheavy corrections using two 
qualitatively-different superpartner spectra, specifically ~ \i = 100 GeV, 
m-g = 300 GeV, m» = 100 GeV, m- q = m i = 1000 GeV and m A = 1000 GeV for the 
"heavy squarks" case, and m~ h ~ /j, = 250 GeV, m g = 300 GeV, rriw = 100 GeV, 
= m £ - = 400 GeV and = 400 GeV for "light squarks." The "cloud" indi- 
cates the region where Ag > 2. These predictions carry estimated uncertainties 
of ~ ±5 GeV from various approximations and from the GUT-scale thresholds 
discussed in the text. Also shown are the estimated allowed mass ranges for the 
running parameter nib as extracted in Sec. V. 
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y G , the splitting of the the superheavy members of the gauge 45 which contributes in 
proportion to a G , and any other superheavy multiplets which may couple to the 10 g 
or the 16 3 . In addition, there is the shift in the prediction of the overall superheavy 
mass scale which occurs if the superpartner mass splittings change the scale at which 
gi and g 2 meet to M' G . This shift could also result from GUT-scale threshold cor- 
rections to gi and g 2 - In any case, such a shift induces an effective GUT threshold 
correction given by Ae„ = £ hx{M' G /M G ) EaiKVc + L a v a G + £ CT K'«y v ) for u = t,b, r. 
After summing over all superheavy particles a, the last term in Ae u becomes equal 
for all v since it only involves complete SO(10) multiplets; therefore we may drop it. 
The first two terms are i hx{M' G /M G ) [(7, 6, 7)y G - (^, 22, ^§)a G ] , which become f3§ 



\n(M G /M G ) (0, —1, 0)2/g — (0, |, after subtracting out an irrelevant constant. 

Then an M' G anywhere between 1.3 x 10 15 GeV and 1.3 x 10 17 GeV changes the top mass 
by less than 2%. 

The threshold corrections proportional to y G arise from the splitting of the 16 3 into 
the right-handed neutrino and the rest of the standard-model matter fields, and of the 
10 H into the superheavy Higgs triplets and the standard-model Higgs doublets. Since 
the right-handed neutrino is but a small part of the 16 3 , its contribution is small, as 
we saw above in Ae„, and since its mass is expected to be at or below the GUT scale 
it can only raise the mass of the top. The splitting of the 10^- is dominated by the 
large hierarchy between the doublets and the triplets; in the language of the Pati-Salam 
subgroup G PS = SU(2) L x SU(2) R x SU(4) of SO(10), we write 10^ -> (2, 2, 1) + (1, 1, 6). 
But in fact any complete multiplet of Gps will contribute equally to \ t , \b and A r , so 
the large splitting between the (2, 2, 1) and the (1, 1, 6) does not generate any threshold 
corrections. Furthermore, the (1, 1,6) decomposes into a 3 + 3 of color SU(3), which as 
it turns out each contribute equally to the three Yukawa couplings (assuming they are 
heavier than the right-handed neutrino), so no threshold contributions result from their 
possible splitting. 

Threshold corrections proportional to a G would be generated by splittings amongst 
the superheavy fields in the gauge 45 as well as amongst the members of the 10 H and 
16 3 . Only splittings of A^ from are not largely suppressed by the fixed-point be- 
havior, as discussed above. The dominant contributions to such splittings come from 
three multiplets 0i,2,3 £ 45, having masses M 1j2i3 respectively, which transform as 
0i ~ (2,3, |), 02 ~ (2, 3, — |) and 3 ~ (1, 3, |) under SU(2) x SU(3) x U(1) Y - They 
yield tf, — e T = {<y G /ir) [ln(M!/M 2 ) — \n(M 3 /M G )]. Even in the somewhat extreme case 
of M 2 ~ 10Mi ~ IOM3 ~ 100M G , the change in m t is only 2 to 6 GeV. 

Finally, there will in general be threshold effects from various other couplings to the 
I63 1Q.H I63 vertex. These obviously depend upon the Higgs content of the theory, and 
cannot be estimated without a concrete model. In particular, if some large multiplet 
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such as a 126 H has a large coupling to the 16 3 and is far from degenerate in mass, then 
large threshold corrections could result; in that case, the Yukawa couplings would not 
in effect be unified. Since we cannot rule out such a possibility (even after we impose 
the restriction that no Landau poles be encountered within, say, and order of magnitude 
of Mq), our predictions will only be valid for models in which either the couplings of 
all multiplets (except the 10 g ) to the 16 3 are small, or such multiplets are practically 
unsplit, or they make equal contributions to Ab and A T . For this last claim we have used 
the fact that \f£ £) (R) + F^(R)\ = \F t (e \R)\ < 1: any threshold effects that are equal 
for Ab and A T or that affect only A t , are greatly suppressed by the fixed-point behavior. 
This applies, for example, to any multiplet that only corrects the 10 H leg of the vertex, 
since it could only distinguish the doublet that couples to down-type fields from the 
doublet that couples to up-type fields but could not distinguish the b from the r. Even 
a 30% threshold correction of this sort would affect the top mass prediction by less than 
3%. 

In summary, while we cannot eliminate the possibility of large GUT-scale threshold 
corrections to our predictions, we have shown that all those corrections which are generic 
to SO(10) SUSY GUTs are not expected to change m t by more than a few GeV. 

X. EXTENSIONS 

The analysis of Sees. II through VIII has actually assumed that the two light Higgs 
doublets lie entirely within a single irreducible representation of SO(10). To what extent 
does the analysis, and the resulting top quark mass predictions, remain valid when the 
Higgs doublets have components in other irreducible representations? In general there 
will be a set of mixing angles {6u,i} and {0o,i\ describing the components of the two 
light Higgs doublets, Hu and H D) in various SO(10) multiplets (10 g , 126 g , or some 
other representation) labeled by i. Suppose for now that the third generation masses 
are generated by a set of Yukawa interactions which are all of the form lf^lO^y 16 3 , 
where the {i'} are a subset of the {i}. In this case the GUT boundary condition Ab = A T 
will occur regardless of the values of {0u,i} and {#D,i}; on the other hand, A t ^ Ab if 
&u,i' 7^ 9d,i' for any i'. The analysis of this paper still applies to this situation provided 
an additional term Ae t is added to e t to reflect this change in boundary condition. As 
we saw above, the fixed-point behavior implies a considerable insensitivity to e t (see the 
corresponding sensitivity function F^). Typically Ae t ~ 0jj^ — 6 D i so for Uti — 9 Dji < 0.3 
the shift in the top mass is less than 5 GeV. (This is not necessarily an uncertainty: in 
a given model the mixings are computable and so is the top mass shift.) This shows 
that our results apply even when several 16 3 10 H 16 3 interactions contribute to the third 
generation masses and Hu and H D have sizeable components in different 10 H multiplets. 
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A particularly interesting subclass of the above models contains just one pair of 10's, of 
which only one, 10 1; couples to the 16 3 . The 10 2 is introduced to make the triplets in 
10 1 heavy This is achieved while keeping the doublets light via the coupling 10 1 45^10 2 , 



provided the A5 H gets a VEV in the B — L direction [0. A more detailed discussion of 
a model of this type will be given elsewhere p(| . 

The restrictions on the values of {#oy} are stricter when i' refers to a 126 H # which 
couples via <^' 16 3 l26 H n 16 3 . because this leads to Ab 7^ A r at the GUT scale. Using our 
analysis then requires an additional contribution to e& — e T of order (£«'/^G)#D,j'- Shifts 
in m t of up to 10 GeV occur when / Xg)Qd,i' ~ 0.1. We conclude that while H v and 
Hp must lie predominantly in a single 10jj, they may have a certain amount of mixing 
with doublets in other multiplets, where the allowed mixing may be quite large if the 
other multiplet is a 10 H but must be small if it is a 126 H with a significant coupling to 
the third generation. 

Another contribution to e^r may arise from operators of the form 16 2 O 16 3 which 
mix the second and third generations. Such operators must be present in order to 
generate V c b, in fact, O must have sufficient structure to give different values for the (2,3) 
entries of the up and down Yukawa matrices. This implies a non-universal contribution 
to £tf>,T- Normally the contribution is of order Uj, ~ 10~ 3 , which hardly affects the top 
quark mass prediction. However, if in some scheme (see, for example, the last reference in 
|i~2fl ) the operator 16 2 O 16 3 is responsible for generating a sizable fraction of m^/m T and 



m s /mt, [rather than having these generated by the (2,2) entries of the Yukawa matrices] 
then the contribution to ej, iT could be large. For example, if there is no operator of the 
form 16 2 O 16 2 to give a mass to the muon or the strange quark, then we would expect 
e b — e T ~ m^/rrir ~ 0.04, which would lead to a shift in m t of ~ 8 GeV. This extreme 
case is however unnatural since such schemes will not lead to an understanding of why 
V c | ~ m c /m t m^/rrir ~ m s /m^ in any case, the shift in such schemes is calculable 
and does not introduce much uncertainty into the prediction of the top mass. 



XI. CONCLUSIONS 

In this paper we have computed the top quark mass in supersymmetric SO(10) grand 
unified theories including all contributions larger than about 5 GeV. The assumption 
underlying this computation is that m t , rrib and m T all originate from a single Yukawa in- 
teraction with an SO(10) multiplet which dominantly contains the two low energy Higgs 
doublets. While this is a somewhat restrictive assumption, it has certain virtues. It 
seems to us to be the simplest assumption, within grand unified theories, which can lead 
to a top mass prediction in terms of low energy quantities. In particular it does not in- 
volve any ansatz about the origin of masses for the lighter two generations. Furthermore 
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this assumption leads to an almost unique group-theoretic structure which underlies the 
top mass prediction: except for a narrow range of (measurable) MSSM parameters, the 
Yukawa interaction must be of the form 16 3 10 H 16 3 . The top mass prediction which 
results from the unification of the three third-generation Yukawa interactions recalls the 
prediction of the weak mixing angle which follows from the unification of the three gauge 
couplings. There is however an important difference between these two cases, having 
to do with threshold effects at the SUSY scale. In the gauge case these effects consist 
of a renormalization of only one operator (the gauge kinetic term), and are therefore 
small as long as the theory is perturbative. In contrast, there are two possible Yukawa 
interaction operators below the SUSY scale, corresponding to the two Higgs doublets, 
that contribute to the mass of each fermion once these doublets acquire VEVs. Thus a 
hierarchy of VEVs could result in large threshold corrections to the lighter masses with- 
out necessarily invalidating perturbation theory. Such a hierarchy is actually mandated 
in the SO (10) unification under our assumptions, and results in another contrast with 
previous unification scenarios, specifically the partial Yukawa unification in SU(5). A 
crucial outcome of the GUT boundary condition Ab = A T is that the resulting value of 
m>b/m T can only be reconciled with experiment if the top quark Yukawa coupling is large. 
This is why in SU(5) and in SO(10) models the top quark is predicted to be heavy. In 
the former, the ratio of A t /Ab at the GUT scale is an arbitrary parameter, which thus 
precludes a prediction of m t . The SO(10) GUT boundary condition fixes A t /Ab to be 
unity; however, we are forced to take vu/vd ^> 1 to account for m t /m bT 3> 1, and 
hence large SUSY threshold effects introduce a new parameter 5mb into the top mass 
prediction. This feature of SO(10) is an improvement over SU(5): now the undetermined 
parameter is a ratio of observable mass parameters, which is measurable in principle by 
future experiments, rather than a ratio of VEVs which will be more difficult to deduce. 
Furthermore, if we assume a hierarchical structure in the SUSY spectrum (a discussion 
of such a structure will be presented elsewhere ||20|| ) in which the sfermions are consid- 
erably heavier than gauginos and higgsinos, there are indeed no large threshold effects 
and the top mass is sharply predicted. 

With these caveats in mind, and choosing the hierarchical spectrum, our basic re- 



sults are shown in Fig. 3 and Table I, and are approximated by Eq. (|32|) and Table II. 
These show the dependence of the predicted top quark pole mass on a 3 and m b and 
the dominant superpartner spectrum effects. Unless future experiments find a value for 
«3 below the range shown in the figure, the hierarchical spectrum predicts a top quark 
heavier than ~ 170 GeV to within 5 GeV or so. (Of course, various uncertainties could 
pile up and result in a lower top mass, but this is unlikely.) A smaller experimental 
value for m t would indicate that Smb is considerable; measuring the various MSSM mass 
parameters could then test this SO (10) unification scenario as outlined at the end of 
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Sec. VI. More generally, our prediction reads r??f ole = A t (177 GeV) [1 + ft{jn>A)\ where 
A t = F(R) + J2 n Fi e) {R)tn + E v F^(R)k v , F(R) and the sensitivity functions F®(R) 
and FW(R) are shown in Fig. 2, the e n and the k u are given in Eqs. ( |23a| - |25cD , and 
R = [ m p(4.1 GeV)/m T ] (Vr/Vb) [1 - a 3 (/iz)/37r] [1 + f R (m A )]. In any specific MSSM- 
based GUT model in which A^ ~ A^ ~ A^, both the possible contributions of the 
diagrams in Fig. 1 and any further deviations from the equality of GUT-scale Yukawa 
couplings can be calculated and inserted into the above expressions. These then yield an 
analytic prediction of the top mass to within a few GeV — which is more than sufficient 
for comparison with top mass measurements likely to be made in the near future. 

Note added: After this manuscript was completed, we received a paper by 
M. Carena, M. Olechowski, S. Pokorski and C.E.M. Wagner |[38|| , in which some of 



the same issues as in our work are addressed and studied in a particular context, namely 
that of universal soft masses at the GUT scale. Many of these issues are discussed in 
detail in our Ref. [F2DJ. A few have also appeared in Ref. |5S . 
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APPENDIX: USEFUL FUNCTIONS 

We collect below various useful functions needed for calculating corrections to m&, R 
and mf ole . The integral which arises in calculating the finite 1-loop corrections to the b 
mass in Fig. 1 is given by 

u du 



r 

I(x,y,z) = / 
Jo 



u + x)(u + y)(u + z) 
xylnx/y + yzlny/z + zx\nz/x 
[{x - y){y - z){z - x)\ 

In some of the SUSY threshold corrections we need the function 
Gi(x,y)= / win [(1 — u) x + uy] du 
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(x-y) 



l — {x 2 (lnz-i) -y 2 (in y- \ 
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2y [x (In x — 1) — y (\ny — . 



(A2) 



When matching R in the 2HSM and the OHSM, we use both I(x,y,z) and Gi(x,y) 
through the function 

IGvr 2 f R (m A ) = (A 2 - A 2 ) Gl o) + \Xl Gl (^f, ^pj 



+ |A 2 Gi 



m A + m z w mt\ 1 2 ( m A + m z w 
, — r 2 A (, i 



/4 



/4 



,0 



+ A t m A I(m A + m W: m t , m w ). 



(A3) 



Finally, for the top quark pole mass, the relevant function is 

/.(>«.<) = + ^ + ¥1 + ¥>* 



where 



and 
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\\ G 2 



Ti. 0±z + a 2 G m * mz 



(1 — w) 2 :r + mj/ d«, 



G 2 (x,y) = [ u 
Jo 

Gs(x,y) = / win — ii) a; — uyl c£u. 
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